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Abstract 

We study integrable geodesic flows on Stiefel varieties V n ,r = SO(n) /SO(n—r) given 
by the Euclidean, normal (standard), Manakov-type, and Einstein metrics. We also 
consider natural generalizations of the Neumann systems on V„ ir with the above metrics 
and proves their integrability in the non-commutative sense by presenting compatible 
Poisson brackets on (T*V n , r )/SO(r). Various reductions of the latter systems are 
described, in particular, the generalized Neumann system on an oriented Grassmannian 
Gn r and on a sphere S n 1 in presence of Yang-Mills fields or a magnetic monopole 
field. 

Apart from the known Lax pair for generalized Neumann systems, an alternative 
(dual) Lax pair is presented, which enables one to formulate a generalization of the 
Chasles theorem relating the trajectories of the systems and common linear spaces 
tangent to confocal quadrics. Additionally, several extensions are considered: the 
generalized Neumann system on the complex Stiefel variety W n , r = U(n)/U(n — r), 
the matrix analogs of the double and coupled Neumann systems. 
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1 Introduction 



A Stiefel variety V nr is the variety of r ordered orthogonal unit vectors (ei, . . . , e r ) in the 
Euclidean space K™ or, equivalently, the set of n x r matrices 

X = (ei---e r ) g M n . r (R) 

satisfying the condition 

X T X = I r , (1.1) 

where I r is an r x r unit matrix (e.g, see In particular, V n \ is a sphere S"™ -1 , while 

both V^, n and V n>n -i are diffeomorphic to SO(n). 

The integrable geodesic flows on V n ,i and V^ lT -, 2 < r < n — 1 are constructed in 
[48] and [10], respectively. The geodesic flows were described in the homogeneous space 
representation SO(n) / SO(n — r) following a general approach to the integrability of geodesic 
flows on homogeneous spaces jTOl HU H2] . 

In the first part of our paper (Sections 2, 3, 4) we study geodesic flows in the redundant 
coordinates X with constraints (jl.ip by using a Dirac approach and the corresponding 
Poisson structure. This allows us to describe the flows in a quite transparent way. Our 
main tool is the theorem on non-commutative integrability of Hamiltonian systems (see 
[391138] ) and for various examples of integrable flows, we calculate the dimension of invariant 
isotropic tori and give the matrix Lax representations (Section 4). 

For the Manakov-type metric on V n ^ r a geometric interpretation of the motion in the 
form of the classical Chasles theorem for the geodesic flow on an ellipsoid is given (Section 
[8]). In addition, the complete integrability of the geodesic flows of the S'0(n)-invariant 
Einstein metrics constructed by Jensen [27] . and Arvanitoyeorgos, Dzhepko and Nikonorov 
[4] are proved in Sections 4 and 6. 

The second part of the paper (Sections 5, 6, 7, 8) considers integrable Neumann flows 
on the Stiefel variety V n ^ r — SO(n)/ SO(n — r) and the oriented Grassmannian variety 
G n .r — SO(n) I SO(r) x SO(n — r), generalizing the famous Neumann system on the sphere 
gn-l_ The latter is defined as a natural mechanical system with the Hamiltonian (see 
[401 133 [36]): 

H neum = -(p,p) + -(Ae,e), A = diag(ai, . . . , a n ), (1.2) 

where the cotangent bundle T*^" -1 is realized as a submanifold of M. 2n {e,p} given by the 
constraints (e, e) = 1, (e,p) — 0. This system, together with the Jacobi problem on the 
geodesic flow on an ellipsoid, provides one of the basic and most beautiful examples of 
application of algebraic geometric tools to integrable systems (e.g, see [Ml 133 HI H] ) ■ It is 
also directly related to many other integrable models (e.g., see [551 l4"9l [501 13"2l |2"2"] ) . 

The Neumann systems on V n ^ r which we consider have the kinetic energy of SO(n)- 
invariant metrics described in Section 3 and the potential function 

V= l -tr{X T AX)= 1 -Y j {e l ,Ae l ). 

i=l 

Two matrix Lax representations are presented. The first, a "big" one, given by Theorem 
15.21 is closely related to the symmetric Clebsch-Perelomov rigid body problem [41] . For 
r = 1, it was given by Moser in [3S] and for r > 1 and the case of the Manakov type 
submersion metrics by Reyman and Semenov-Tian-Shanski [44] within the framework of 
the i?-matrix method. Note that for r > 1 this Lax pair does not define a Neumann system 
on V n>r uniquely and does not provide a non-commutative set of integrals, necessary for the 
integrability. 
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In contrast, the second (dual, or "small") Lax pairs, given by Theorem 18 . 1 1 are equivalent 
to the Neumann systems with the Euclidean and normal metrics up to an action of a finite 
discrete group. For the Neumann system with the Euclidean metric, the small Lax pair was 
first given in the unpublished manuscript [28] . 

In Sections 6 and 7 we give a detailed proof of the non-commutative integrability of 
the considered Neumann flows by using the Bolsinov completeness condition for a set of 
Casimir functions of the pencil of compatible Poisson brackets (see [9]). We also indicate an 
integrable generalization of the Neumann system on Grassmannians with a quartic potential. 

In Section [5] we propose a geometric interpretation of the integrals of the Neumann 
systems on V n>r obtained from the dual Lax representation. Our geometric model generalizes 
the celebrated Chasles theorem adopted by Moser for the case r = 1 (see, e.g, Theorem 4.10 
in [33). 

Magnetic Neumann flows in the rank two case (on G n ^2 and V n ^) as well as the motion 
of a particle on a sphere S 71 ^ 1 under the influence of a Yang-Mills field are presented in 
Appendix 1. Finally in Appendix 2, we briefly consider the rank r double and coupled 
Neumann flows, as well as an extension of the Neumann system onto a complex Stiefel 
variety W n , r = U(n)/U(n — r). 

The geodesic flows and Neumann systems considered in this paper are written in a form 
appropriate for their integrable discretizations, which we describe in a separate paper [23j . 

2 Hamiltonian Flows on Stiefel Varieties 

Stiefel varieties. As it was mentioned in Introduction, a Stiefel variety V n , r is the vari- 
ety of r ordered orthogonal unit vectors (ej., . . . , e r ) in the Euclidean space (R™, (•, •)), or, 
equivalently, the set of n x r matrices satisfying constraints (jl.ip . Thus V n , r is a smooth sub- 
variety of dimension N = rn — r(r + l)/2 in the space of n X r real matrices M n ^ r (M.) = M. nr 
and the components of X arc redundant coordinates on it. 

The left SO(n) action on V n ^ r (X i-> RX, R € SO(n)) is transitive, hence V n , T can be 
realized as a homogeneous space of the Lie group SO(n) as well. If fix the orthonormal base 
in R" 

Ei = (1,0,,0,...,0) T , E 2 = (0,1,0,... ,0) T , E n = (0,0,0,. . . , 1) T (2.1) 

and take the point X = (E\, . . . , E r ) € V n<r , then the orthogonal transformation fixing Xq 
(relative to the above basis of R") must have the form 

^ BeSO(n-r). (2.2) 

Since the isotropy group of Xq is isomorphic to SO(n — r), the variety V n , r can be identified 
with SO(n)/SO(n-r). 

The Poisson structure. The tangent bundle TV n ^ r is the set of pairs (X, X) subject to 
the constraints 

X T X = I r , X T X + X T X = 0. (2.3) 

On the other hand, the cotangent bundle T*V n , r can be realized as the set of pairs of n x r 
matrices (X, P) that satisfy the constraints 

X T X = I r: X T P + P T X = 0. (2.4) 

The latter give r(r + 1) independent scalar constraints 

Fij = (ei, e,) - 5 tJ = 0, = (ei,pj) + (e jl p l ) = 0, 1 < i < j < r, (2.5) 
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where pj is the j-th column of the matrix P. This realization of T*V n , T is motivated by 
description of the geodesic flows of the Euclidean and normal metric on V ntT given below, 
however there are other natural realizations of T*V n ^ r (see Section 13]). 

The canonical symplectic structure u> on T*V n . r is the restriction of the canonical 2-form 
in the ambient space T*M„ jr (]R), 



U) 



dfi df 2 \ 






dei ' dpi J 




de t 



For our purposes it is convenient to work with the redundant variables (X, P) . The 
canonical Poisson structure on T*V n ^ r can then be described by using the Dirac construction 
[3J[Tni[3S]. Namely, let {•, -}o be the canonical Poisson bracket on M. 2nr 

r 

{/i,/ 2 }o = E 

and Cij be the inverse of the matrix {i^, Fj}o, i, j = 1, . . . , r(r + 1), where, for the sake 
of simplicity, we denoted constraints (|2.5p by Fi = 0, i = 1, . . . , r(r + 1). Then the Dirac 
bracket is given by 

{h, /a} - {/i,Mo + Y^{r,.I:\j;,{i;-hU. (2.6) 

The subvariety T*V n . r appears as a symplectic leaf of the Dirac bracket and the restriction 
of {/i, f^} to T*V ll:r depends only on the restriction of /i and to T*V n . r . 
The Hamiltonian equation 

f = {f,H} 

can be also described by the using Lagrange multipliers. We shall write them in the matrix 
form: 

p = ~ ox- + XA + pu > 

where A and II are r x r symmetric matrix Lagrange multipliers uniquely determined from 
the condition for the trajectory (X(t), P(t)) to satisfy constraints (|2.4[) . 



Momentum mappings. The Lie group SO(n) naturally acts on T*V n ^ r by left multipli- 
cation: 

R ■ (X, P) = (RX, RP), ReSO(n). (2.8) 

Below we use the well known identification of A 2 M™ with a subset of so(n): x A y — 
x <S> y — y <S> x = x ■ y T — y ■ x T , x, y £ K™. Also, (•, •) is proportional to the Killing metric 
on so(n): 

= -^tr(&6). (2-9) 

^1:^2 £ so(n). By the use of the scalar product (|2 .9[) we identify so(n) and so(n)* . 

Proposition 2.1 The left SO(n)-action \2. 8\) is Hamiltonian. The equivariant momentum 
mapping $ : T*V n _ r — > so(n)* = so(n) is given by 

r 

$(X,P)=J2Pi Ae i5 (2.10) 

i=l 

or, in the matrix form $(X, P) = PX T — XP T . 
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Proof. The left S0(n)-action (|2~%|) on (T*M n>r (R), {•, -} ) is Hamiltonian with the momen- 
tum map (|2.10[) . The constraint functions X T X and X T P + P T X are S'0(n)-invariant. 
Therefore 

{^,^10 = 0, i = l,...,r(r + l), (2.11) 

where $^(X, P) = P),£) is the Hamiltonian function of the action of the one- 

parameter subgroup {exp(s£),s G M}. 

In view of the definition of the Dirac bracket (|2 .6[) and (|2 . 1 1|) . the Hamiltonian flows 
of on T*V nyr with respect to the brackets {■, -}q and {•,•} coincide. This proves the 
proposition. □ 

Together with a left 5'0(n)-action, we also have the natural right S'0(r)-action: 

(X,P)-Q = (XQ,PQ), QeSO(r). (2.12) 

Following similar lines, one can prove 

Proposition 2.2 The right SO(r)-action H2.12\) is Hamiltonian. The equivariant momen- 
tum mapping ^ : T*V n . r — > so(r)* = so(r) is given by 

m{X 1 P) = X T P-P T X. (2.13) 

The momentum mappings $ and 'J are Poisson with respect to the (+) and (-) Lie- 
Poisson brackets on so(n) and so(r): 

{hi o P), h 2 o *{X, P)} = h 2 (ji)}+ ln) , t l = *(X, P), 

{h ° P), h ° P)} = ih(v), f2(v)} Mr) , v = *(x, P), 

where 

{ht(n), h2(ji)}+ (n) - (M, [V/i x (m), V/i 2 (m)]), ^1. ^2 : so(n) -> R, 
{/iW,/2(r?)}- (r) = -ft, [V/ifa), V/afa)]), /i,/ 2 : «o(r) -> K 
and where the brackets (•, •) denote the scalar product (|2.9[) on so(n) and so(r), respectively. 

The algebra of 50(n)-invariant functions. Let C°° (T*V n . r ) so ^ be the algebra of 
50(n)-invariant functions on T*V n , r . Since SO(n) acts in a Hamiltonian way on T*V n>r , 
C°°(T*V n!r ) so ^ is closed under the Poisson bracket. 

Let Xq = (Ei, . . . , E r ). The S'0(n)-invariant functions, via restrictions, are in one-to-one 
correspondence with the SO(n — r)-invariant functions on T£ V n , r - 

We can write Po G T^- g V n ,r as a block matrix 

Pa = (g) . (2-14) 

where Pi and Pi are r x r and (n-r)xr matrixes, respectively. Then Po satisfies constraints 
(pO|) at X = Jf if Pf = -Pi. Also, the SO(n - reaction on TJ o K, r is given by 

Px>—>Pl, P 2 ^B-P 2 , BeSO(n-r). (2.15) 

Lemma 2.3 The maximal number of functionally independent SO(n) -invariant functions, 
i.e., the differential dimension of C°° (T*V n-r ) so ^ n \ equals 



MimC°°(T*V n , r ) so{n) 



dimV^.r — dimS'0(ri — r), n <2r 

dim V n r — dim SO(n — r) + dim SO(n — 2r), n > 2r. 

(2.16) 
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Proof. The differential dimension of C°° {T*V rhr ) sotKn ^ is just the codimension of the generic 
SO(n — r) orbit in V n , r . The dimension of the orbit SO(n — r) ■ Pq is 

dim SO(n — r) — dim SO(n — r)p , 

where 5(3(n — r)p is the isotropy group of Pq. Since SO(n — r) Po = {I„_ r } for n < 2r, 
for a generic Pq we get the first relation in (I2.16[) . Further, in the case n > 2r, we can take 
P = (E r+u E r+2 ,...,E 2r ). Then SO(n-r) Po = SO(n-2r). □ 

It can easily be verified that the restrictions of the S'0(n)-invariant functions 

= {euPj) - (ej,Pi), (P T P)ij = (Pi,pj), i,j = !,•••, r (2.17) 

to T^ o V ntr define the generic orbits of the action (|2.15[) . In particular, we get the following 
simple statement. 

Lemma 2.4 If a smooth function f G C°°(T*V n ^ r ) Poisson commutes with functions 
then it Poisson commutes with all SO(n) -invariant functions on T*V n>r . 

The Poisson bracket on C°° (T*V n ^ r ) so ^ can be described as follows. The restriction 
of momentum mapping $ to T^ o V n . r establish the isomorphism 

Tx Vn,r = », (2.18) 

where is the orthogonal complement of so(n — r) in so(n). Within identification (I2.18p . the 
SO(n — r)-action (|2.15[) corresponds to the adjoint SO(n — r)-action on and the Poisson 
bracket on C 00 (T*V n _ r ) so ^ corresponds to the Poisson bracket 

{A(0,/2(0}» = [v/i(e), V/ a (0]). 

on the algebra C ,00 (o)' SO( -™~ r ) of SO(n — r)-adjoint invariants on U (see Thimm [48]). 

Recall that the set of commuting 5 l O(n)-invariant functions 21 is complete if it contains 
maximal possible number of independent functions, that is (see [TUl IT2"] ). 

ddimSl = &va\V ntr — I, 

where 21 is the dimension of a generic adjoint orbit in 

$(T*V n , r ) = Ad SO („) $(?x V n ,r). 

Note that, for n < 2r + 1, the generic orbit in &(T*V, lyr ) is regular, while it is singular 
otherwise. One can prove the following relations. 

Lemma 2.5 

2/= / ^-[f], n<2r + l, 
\ 2r(n-r-l), n>2r + l. 

Now consider the chain of subalgebras so(n — r + 1) C so(n — r + 2) C • • • C so(n), where 
a matrix £ £ so(n — r + i) is included in so{n) as a block matrix 




£ 



(2.19) 



Let 2ti be the algebra of invariants on so(n — r + i) considered as a polynomials on so(n) 
and restricted to 0. Then 

2t = 2li + ---+2t r (2.20) 

is a complete polynomial commutative subset of C' 00 (t)) ,so ( n-r ) (see again |T0j [12] ) ■ Other 
complete commutative sets of SO (n)-invariant functions are given in |17j and Theorem 16.41 
below. 
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3 Geodesic Flows 



S'0(n)-invariant metrics. An 5 l O(n)-invariant metric g on V n ,r can be specified by a 
positive definite, SO(n ~ r)-invariant scalar products g$ at the point X as follows: 

g(R-r)x,R-r)2)x = go(r)i,V2), V2 G T Xo V nir , X = R-X . 

Equivalently, an 5'0(n)-invariant metric can be denned by using an S l O(n)-invariant 
Hamiltonian function H(X,P), which is quadratic in momenta and positive definite on 
T*V 

There are two natural S l O(n)-invariant metrics on the Stiefel variety V n y. the metric 
induced from the embedding of V n , r in the Euclidean space M„ )r (K) and the normal metric 
induced from a bi-invariant metric on the Lie group SO(n) (see below). 

Concerning geometrical significance, one should also mention invariant Einstein metrics 
(see [5]) constructed in [27], [I]. 

The Euclidean metric. The Euclidean metric in M njT .(R) is given by the Lagrangian 
function 



L E (X,X) = \tT{X T X) = \^{e u < 
The Legendre transformation 



2 

t=i 



p= dLE =jc (31) 
dX 



yields the Hamiltonian function 

H E (X,P) = \tT(P T P) =J2l(Pi,Pi) (3-2) 

i=l 

defined on the cotangent bundle T* M n ^ r (M.). 

We shall refer to the restriction of the above metric to V n ,r as the Euclidean metric, which 
will be denoted by ds%. The geodesic flow is described by the Euler-Lagrange equations 
with multipliers 

d L d_LE = dL R - = k = _ ± t ± 

dt dX dX ' 

where the symmetric matrix A is uniquely determined from the condition for the trajectory 
X(t) to satisfy the constraints X T X = I r . 

Taking into account constraints (|2.3I) and the Legendre transformation (|3 . 1 1) , we see that 
the cotangent bundle T*V n ^ r can be represent as a submanifold of T*M n r (R) given by (12.41) . 
The corresponding Hamiltonian flow of He(X,P) with respect to the Dirac bracket is 

X = P, 

T (3-3) 
P = -XP T P. 



Submersion metrics. Let g(-,-) be a right-invariant metric on SO(n). The subgroup 
SO(n — r) acts freely on SO(n) by isometries (right action). There is a g-orthogonal de- 
composition of TjiSO(n) 

SO{n-r) B . + n R = T R SO{n), ReSO{n), 
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where SO(n — t)r is the tangent space to the fibre R ■ SO(n — r). By definition, the 
submersion metric <?(-, ■) on SO(n)/ SO(n — r) is given by 

ff(fx,6)p(ii) = ff(6b6)fl, 6 e T p(R) (SO{n)/SO(n - r)), & G 6 = 

where p : SO(n) — ► SO(n)/SO(n — r) is the canonical projection (e.g., see [S]). The vectors 
in 50 (rt — r)^ and Hr are called vertical and horizontal respectively. 

The Hamiltonian of a right invariant metric on SO(n) can be written in the form lio$, 
where h is a positive definite quadratic form on so(n) and $ : T*SO(n) — > so(n) is the 
momentum mapping of the natural left SO (n)-action. It follows that the class of submersion 
metrics on SO(n) / SO{n — r) is given by Hamiltonian functions of the form h o $, where 
now $ is the momentum mapping of a natural left S'0(n)-action on T*(SO(n)/SO(n — r)) 
(e.g., see [14]). 

The above observation helps us to write down the Hamiltonians and the geodesic flows 
of the submersion metrics in the representation of the cotangent bundle T*V n>r given by the 
constraints (|2.4[) . The Hamiltonians are 

H A (X,P) = h A o $ = i(A$, $) = ~tr (A(PX T - XP T )(PX T - XP T )) , (3.4) 

where ft&(£) = £)> £ e so(n) and A : so(n) — > so(n) are positive definite operators. 

Proposition 3.1 The equations of the submersion metrics geodesic flow generated by {3.$ 
are 

X = A($)-X, 

(3.5) 

P = A($) ■ P. 

Proof. By using the chain rule dH = dh A o c?<i>, one gets expressions (|3.5p for the derivatives 
of the Hamiltonian (|3.4[) with respect to P and —X. Further, it can easily be verified that 
then the Lagrange matrix multipliers A and n in (12.71) are zero. □ 

Remark 3.1 In particular, for r = n, the system (|3.5I) describes the right-invariant geodesic 
flow on the Lie group SO(n). The symmetric form of the equations differs from the sym- 
metric representation of the rigid body equations given in [6]. 

The normal metric. If g(-, •) is a bi-invariant metric induced by the scalar product (|2.9[) . 
then the submersion metric is called the normal metric. It is proportional to the standard 
metric induced by the negative Killing form 0) on SO(n)/SO(n — r). We shall denote the 
normal metric by dsg. Contrary to a generic submersion metric, the normal metric is also 
S'0(n)-invariant. The corresponding Hamiltonian has the form 

H (X,P) = = \tr{P T P) Ur((X T Pf), (3.6) 

and its geodesic flow is given by 

X = $ • X = P - XP T X, (3.7) 
P = $ • P = -XP T P + PX T P. (3.8) 



Under the conditions (|2.4I) . the relation (|3 . T[) can be uniquely inverted and one gets 

P = X - ^XX T X. (3.9) 
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Relations (I3.7[) . (|3.9|) give identification of TV n , r and T*V n . r by means of the normal metric. 
Therefore, the Lagrangian function for the metric c?s§ is 

L Q (X,X) = ^ a ^ ), (3.10) 

where $ (X,X) = $(X,P(X,X)): 

$o(X, X) = XX T - XX T - ^X[X T X - X T X]X T . (3.11) 

For a right S f O(?i)-action (and therefore with opposite signs in the equations), the relation 
(|3.1ip is established in [25] by studying nonholonomic LR systems on the Lie group SO(n). 
Note that the realization of T*V n ^ r via (|2.4p is also natural if we consider the flow given by the 
normal metric. Namely, defining momenta by the Legendre transformation P — OLq/OX, 
from the constraints (|2.3|) we get the conditions (|2.4p . 

Einstein metrics. The momentum mappings $ and \& are invariant under the SO(r) and 
SO(n) actions, respectively. Therefore, the Hamiltonians of the form 

H K (X,P) = + = \tx{P T P) Q + tr((X r P) 2 ) (3.12) 

is SO(n) x 50(r) invariant. Within the class of the metrics ds 2 . determined by the Hamil- 
tonian functions (|3.12[) there is the normal metric (k = 0) and the Euclidean metric 
(k = —1/2). Moreover, for r = 2 there is a unique value of k, while for r > 2 there 
are exactly two values, such that dsf. is an Einstein metric (see [27l H]). Following [4], we 
refer to these metrics as the Jensen metrics. 

Further, in [4], new examples of the Einstein metrics are obtained within the class of 
metrics that we shall describe below. 

Consider the Lie subalgebra 

so(ri) © so(r 2 ) © • • ■ © so(r k ) C so(r), r\ + r 2 H h r fc = r, 

naturally embedded into so(r). Define the Hamiltonian Hk via 

H K = -($,$) +i<To*, (3.13) 

where the quadratic function X is 

* (0 = + y<Ci,a> + • • • + f (3.14) 

Here & are orthogonal projections (with respect to the so(r)-Killing metric) to so(r,) C 
so(r). 

In [1] Arvanitoyeorgos, Dzhepko, and Nikonorov proved that if r% = r% = ■ ■ ■ = r^, 
«i = K2 = • • • = tfc, k > 1, n — r > r\ > 3, than among the metrics defined by Hamiltonians 
()3.13p there are four Einstein metrics (two of them, with k\ = 0, are the Jensen metrics). 

S l O(r)-invariant metrics. Let A be a symmetric, positive definite n x n matrix. The 
geodesic flows on V n , r with Lagrangians of the form 

L A (X,X) - Ur{X T AX) = ^Y / (e i ,Ae i ), 
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were studied in [7J from the point of view of the variational and optimal control problems. 
The Lagrangian La can be considered on the whole space M nxr (R), where the Legendre 
transformation 

P= ^A = AX, (3.15) 
dX 

gives the following Hamiltonian on the cotangent bundle T* M„ )r (K) 
H A (X,P) = ltr(P T A- 1 P) = i^^-Vi). 

»=1 

^From (|2.3[) and (|3.15p we conclude that the cotangent bundle T*V n ^ r can be represented 
as a submanifold of T*M n>r (M.) given by the equations 

X T X = l r , X T A- 1 P + P T A~ l X = 0. (3.16) 

Again, one defines the Dirac bracket with respect to the constraints (I3.16[) . 

Let ds\ be the metric defined by the Lagrangian La and {•, -}a be the new Dirac bracket. 
Then the geodesic flow of ds A can be described by the Hamilton equations 

f = {f,H A (X,P)} A 

restricted to the symplectic leaf Q3.16p . The corresponding matrix form of the flow is 

X = A~ l P, 

(3.17) 

P = XA, 

where A is a symmetric r x r matrix uniquely determined from the condition 

AX T A~ 1 X + X T A~ 1 XA + 2P T A~ 2 P = 0. 

Note that for r = 1, the metric ds 2 A is a standard metric on the ellipsoid (x, A~ x x) — 1, 
while for r it is a right-invariant Manakov rigid body metric on SO(n) (see 35 ). Further- 
more, ds\ is right S'0(r)-invariant and, via submersion, induces a metric on the oriented 
Grassmannian variety G„ ir (see Section [Jj). It would be interesting to prove integrability of 
the corresponding geodesic flows. 



4 Integrability of Geodesic Flows 

The normal metric. As shown in [TQl [IT], the geodesic flows of the normal metrics ds§ 
on the homogeneous spaces G/H of compact Lie groups G are completely integrable in the 
non-commutative sense. The proof is based on the following geometrical observation. Let 
(M, G, $) be a Hamiltonian G-space with an equivariant momentum mapping $ : M — > g*, 
where G is a compact group. Consider the following two natural sets of functions on M: the 
functions obtained by pulling-back the algebra C°°(g*) by the moment map and the set of 
G-invariant functions C°°(M) G . They are closed under the Poisson bracket and according to 
the Noether theorem {$*(G°°(g*)), G°°(M) G } = [26]. Moreover, $*(G°°(0*)) + G°°(M) G 
is a complete set of functions [IT]. That is, any Hamiltonian system with those integrals 
is non-commutatively integrable [321 13H]- In particular, consider the case when (M, G, $) 
is a cotangent bundle T*{G/H) with a natural G action. Since the Hamiltonian Hq of the 
normal metric is of the form h o $, where h is an invariant quadratic polynomial on g*, 
the function H Poisson commutes with all G-invariant functions (the Noether theorem), 
as well as with all the functions in $*(G°°(g*)) (the mapping h M> h o $ is a morphism of 
Poisson structures). Therefore, the flow of H is non-commutatively integrable. 
Let 21 be the dimension of a generic orbit in <&(T*V ntr ) (see Lemma \2.h\i . 
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Theorem 4.1 [TU] [TT] The geodesic flow of the normal metric j3.7\ ), 13.8)) is completely 
integrable in the non- commutative sense. The complete algebra of first integrals is 

<t>*(C°°(so(n)) + C°°{T*V n , r ) so[n) . 

The generic motions of the system are quasi-periodic over the isotropic tori of dimension 

S Q = 2dimK,r - ddimC , °°(T*y n , r ) so( " ) - 21. 

Recall that the number of functionally independent S'0(n)-invariant functions is given 
by png]) . Thus we have 



So 



r , n > 2r 

dim V„, r + dimS'0(n — r) — dimS'O(rj) + rank 50(7?,) , n < 2r 



It is interesting that for n > 2r the dimension of the invariant tori coincides with 
the dimension of invariant tori of geodesic flows of normal metrics on the corresponding 
Grassmannian manifolds G n r . 



Manakov metrics. The above construction has a natural generalization to a class of 
geodesic flows of submersion metrics given by the Hamiltonians (|3.4[) . such that the corre- 
sponding Euler equations on so(n): 

f = UM} + so(n) £=[V/ia(0,£] = [A&£], £eso(n) (4.1) 

are completely integrable. For example, choose the Manakov operator 

A(Ei A Ej) = bl - hj Ei A Ej, 1 < i < j < n, i.e., A(£) = ad^ 1 ad B ^, (4.2) 

Cb'l Qjj 

where all the eigenvalues of A = diag(ax, . . . , a„) and B — diag(6 l7 . . . , b n ) are distinct and 
A is positive definite. Then the Euler equations (|4.ip are completely integrable [2D] ). 
Moreover, for generic A, among the integrals tr(£ + XA) k one can always find a complete set 
of commuting integrals hi,,.., hi on a generic adjoint orbit in $(T*V^ >r ) (see Brailov [T4"] 
and Bolsinov [9]). Therefore, according to Theorem 2.2 in [TT], we have 

Theorem 4.2 The geodesic flow of the submersion metric 113. 5\) with the metric given by 
the Manakov operator J^.i?| ) is completely integrable in the non- commutative sense with a 
complete set of polynomial integrals \2.17^ and 

tr ^XA + Y^Pi Aejj , k=l,...,n, A e R. (4.3) 

The Generic motion of the system is quasi-periodic over the isotropic tori of dimension 

S = 2dim V„,r - ddimC°°{T*V n:r ) so{n) - I. 

As in the case of geodesic flows of normal metrics, for n > 2r the dimension of generic 
invariant tori is simply S — r(n — r). 

Remark 4.1 The non-commutative integrability implies the usual commutative integrabil- 
ity, at least by means of smooth commuting integrals [TT] . For the case of the geodesic flows 
considered above, the commuting integrals can be taken to be the polynomials (|2.20j) and 
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Like in the case of the right-invariant metric on SO(n), the Manakov metric on V n ^ r 
possess a matrix Lax pair. 

Theorem 4.3 Equations i3. 5\) with the metric given by the Manakov operator imply 
the matrix equation with a spectral parameter A 

C man (X) = $ + AA, A man {\) = A($) + \B . 

The dual Lax pair for the Manakov flows. Consider the Manakov operator (|4.2I) for 
B = A 2 . Then 

A(0=A£ + £A, (4.4) 

and equations (|3.5j) become 

X = A(PX T - XP T )X + (PX T - XP T )AX, 

P = A(PX T - XP T )P + (PX T - XP T )AP. ' ' " 

Theorem 4.4 Up to the action of a discrete group generated by reflections with respect 
to the coordinate hyperplanes in 1", 

(X,P)^(SiX,SiP), i = l,...,n, (4.6) 
Si(xi,...,x n ) = (yi,...,y n ), Vj=xj, j^i, Vi = -xi, 

the geodesic flow |^.5| ] is equivalent to the matrix equation 

^^mmi(^) = [^manWi^mant^)] (4-7) 

with a spectral parameter A and 2r x 2r matrices 

,* iy\ _ (—X T (\l n - A)~ 1 P -X T (\I n ~A)-iX\ 

*-m«nW - ^ pT _ jy-lp pT _ A) -l X J , (.4.8) 

a* (x] _ ( X T (A + \l n )^P X T (A + XI n )^X\ 

•AmanW ~ I _pT (A + XI^P -P T (A + XI^X ) ' W 



Note that after imposing the condition X T P = 0, equations (|4.5j) formally coincide with 
the equations describing rank r solutions of the Manakov system on so(n) (see [2Tj). 

The Euclidean and Jensen's metrics. Since the Hamiltonian Q3.12p is SO(n) x SO(r)- 
invariant, we can apply the general construction used in Theorem 14.11 with respect to the 
SO(n) x 50(r)-action (see [TT]). Let c^(T*K,r) 50( " )><S0(r) be the algebra of SO{n) x 
S'0(r)-invariant functions on T*V n , r - 

Theorem 4.5 The geodesic flows of Jensen's metrics ds 2 . with the Hamiltonian functions 
i3.12\) are completely integrable in the non-commutative sense. The complete algebra of first 
integrals is 

$*(C°°(so(n)) + **(C°°(so(r)) + C°° (ri4,r) SO(re)xSO(r) . (4.10) 
In particular, the geodesic flow \3.3\) of the Euclidean metric is completely integrable. 

The complete commutative set of polynomials # within C°° (T*V rlyr ) S0( - n ' >xS0 ^ as well 
as the integrability of the geodesic flows with Hamiltonians (|3.13p will be given below (see 
(1O0)) in Section 6, Theorem IQ1 and Corollary loT5l) . 
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Remark 4.2 Both the geodesic flows of the Euclidean and the normal metric share the 
isotropic foliation defined by integrals (|4.10[) , but do not share the isotropic foliation defined 
in Theorem 14.11 Namely, the straightforward calculations show that functions (j£>i,Pj) in 
(|2.17|) do not Poisson commute with (iff, iff), and, therefore, the algebra of S'0(n)-invariant 
functions is not conserved along the geodesic flow of the Euclidean metric (see Lemma [2~4]) . 

Here we note the following characterization of SO(n) x SO(r )-invariant metrics. 

Proposition 4.6 If the metric ds 2 on V n . r is SO(n) x SO(r) -invariant, then, up to multi- 
plication by a constant, it coincides with ds 2 . for some k. 

Proof. The statement follows from the fact that the restriction of the Hamiltonian function 
to T^ g V n ^ r is a quadratic form invariant with respect to the transformations (|2.15p and 

Pi i— ► Q^PiQ, P2 i ^ P2Q, Q G SO(r), 
where Pi and P2 are defined in (|2. 14[) . □ 



5 The Neumann Systems on Stiefel Varieties 

The celebrated Neumann system on a sphere S 1 ™ -1 is defined as a natural mechanical system 
with the quadratic Hamiltonian f| 1 . 2[) . By analogy, we define a Neumann on the Stiefel 
variety V n>r as a natural mechanical system with an SO(n) -invariant kinetic energy and the 
quadratic potential function 

V= l -t?{X T AX)= l -Y.{e ll Ae l ). (5.1) 

i=l 

Note that the above potential is constant for r. 

While on the sphere S*"" 1 an S'0(n)-invariant kinetic energy is unique (up to multi- 
plication by a constant factor), on the variety V n ^ r with r > 1 there are many different 
50(n)-invariant metrics. Following Sections 3 and 4, we consider the kinetic energy deter- 
mined by the metrics ds 2 K (see eq. (|3.12p ). Thus, the Hamiltonian has the form 

H neum , K (X, P) = ~ tr(P T P) - ( 1 + k) ti{{X T P) 2 ) + ~ ti{X T AX) . (5.2) 



(5.3) 



2 v V 2 / 2 

Proposition 5.1 The Neumann system with Hamiltonian 15.2(1 is given by 

X = P - (1 + 2k)XP t X, 

P = -AX - XP T P + (1 + 2k)PX t P + XX 1 AX. 
Proof. It is a direct calculation. We have 
f) J-f 

x = - — xn = p - (1 + 2k)xp t x - xii, 

oP 

f) J-f 

P = -— + XA + PU = -AX + (1 + 2k)PX t P + XA + PU. 
uX 

Differentiating the constraints (|2.4|) with respect to time gives 

X T X + X T X = 0. X T P + X T P + P T X + P T X = 0. 
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The first equation implies that the Lagrange multiplier LT equals zero, while the second one 
yields A = X T AX — P T P. ' ' □ 

Note that Hamiltonians (|5.2I) are right 5 l O(r)-invariant, so the momentum mapping 
(|2.13p is conserved by the flows (|5.3[) for any parameter k. In particular, for k = we get 

the Neumann system with the normal metric given by 

X = P - XP T X, 

P = -AX + PX T P + XA = -AX + PX T P - XP T P + XX T AX, ' 

while for k = —1/2 we get the Neumann system with the Euclidean metric with the corre- 
sponding Hamilton equations 

X = P, 

(5 5) 

P = -AX + XA = -AX - XP T P + XX T AX. 

The Lax pair. Although for different n the flows (|5.3[) do not coincide, the derivatives of 
the momentum $ and of the symmetric matrix XX T are the same: 

±* = [XX T ,A], j t (XX T ) = ^,XX T }. (5.6) 
As a result, the following theorem holds. 



Theorem 5.2 Equations \5.3\) . in particular {5.$ and \5.5\) , imply the same n x n matrix 
Lax representation with a spectral parameter X: 

Lneum{\) = A$ + X X T - \ 2 A, A neum {\) = $ - XA. (5.8) 
The proof is immediate. The coefficients of the spectral curve 

T C C 2 {A, v) : det(£„ eum (A) - /il„) = (5.9) 
give us the commuting integrals of both systems, which can be expressed in the form 

3 = {tr{X{PX T - XP T ) + XX T - X 2 A) k | k = 1, . . . , n, X G R}. (5.10) 



The Lax representation (|5.7[) is closely related to the Clebsch-Perelomov rigid body 
system [41]. For r = 1 it was given by Moser in [35] and for r > 1 in [44], without giving 
explicitly the equations of motion. (As was mentioned above, the Lax pair does not define 
the system itself.) The book [44] also describes the Neumann flows on Grassmannians G n . r , 
as well as the magnetic Neumann flow on G n ,2- We shall consider these systems together 
with the magnetic Neumann flows on V n ^ in detail in Section 7 and Appendix 2, respectively. 

Alternative (or dual) Lax pairs, which are does equivalent to equations (|5.5p . (I5.4p (up 
to the action of a finite discrete group) are given below in Section |U 

6 Compatible Poisson Brackets and Integrability 



The Neumann systems on V n ^ r admitting Lax pairs with the Lax matrix (|5.8[) can be ob- 
tained as appropriate reductions of integrable n-dimensional tops having symmetric inertia 
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tensors and moving in a force field with a quadratic potential (the Bogoyavlenski general- 
ization of the Clebsch-Perelomov system [8]), see Appendix 1 and However, these are 
reductions on singular coadjoint orbits in the corresponding Lie coalgebras, hence the inte- 
grability of the Neumann systems does not follow directly from that of the n-dimensional 
tops. 

We now prove a non-commutative integrability of systems (|5.3j) by using the Bolsinov 
condition for a set of Casimir functions of the pencil of compatible Poisson brackets to be 
complete [9]. 



Theorem 6.1 Let all the eigenvalues of A be distinct. Then the Neumann systems \5. 3\) . 
in particular {5.$ and L5.5\) . are completely integrable in the non- commutative sense with 
the non- commutative set of integrals given by \5.1U\) and by the components of the SO(r)- 
momentum mapping h2.13\) . The generic trajectory (X(t), P(t)) corresponding to the maxi- 
mal rank of the momentum is quasi-periodic over isotropic tori of dimension 

if, , r(r — 1) [r~\\ rr~\ 

- 2 {^n-r) + -L_J -[-]) + [-]. (6.1) 

Proof. First, we give the interpretation of the integrals (|5. 10[) from the bi-Hamiltonian 
point of view. Consider the Lie algebra gl{n) of n x n real matrices equipped with the scalar 
product (X, Y) = — htr{XY) and the orthogonal decomposition gl(n) = so(n) + Sym(n) 
onto skew-symmetric and symmetric matrices: 

[so(n), Sym(n)} C Sym{n), [Sym(n), Sym{n)] c so(n). (6-2) 

The scalar product (•, •) is positive definite on so(n) and negative definite on Sym(n). 

Let us identify gl(n)* and gl(n) by means of (•, •). On gl(n) we have a pair of compatible 
Poisson brackets given by the following Poisson tensors 

Ai(€ + 77, C + = (*, K. C] + K, 0} + fa C]>, 

Mt + v,( + 0)\ x = {x-Aiz+r),<; + o}), ' 

where x G gl(n), C G so(n), r},6 G Sym(n) (see 43, 9]). The tensor Ai corresponds to the 
canonical Lie-Poisson brackets in the dual to the semi-direct product so(n) © a d Sym(n). 
Consider the Poisson pencil 

A AljA2 = AiAi + A 2 A 2 , n = {A Al , A2 | Ai, A 2 e M, \\ + \\ ± 0}. 

For A 2 7^ 0, the bracket A All A 2 is isomorphic to the canonical Lie-Poisson bracket on gl(n) 
and its Casimir functions have the form 



where x — h + v, h £ so(n), v G Sym(n) (see (43j [9])- Let 

T = {tr(Aft + « - A 2 A) fc | k = 1, 2, . . . , n, A G M}. (6.4) 

be the union of all the Casimir functions of the brackets with Ai + A 2 ^ 0, A 2 ^ 0. Then F 
is a commutative set with respect to all the brackets in H and, if the eigenvalues of A are 
distinct, J 7 is a complete commutative set on a generic symplectic leaf in (gl(n), Ai) (see 
Theorem 1.5 in [5]). The mapping 

r r 

6 = $ + XX T = ^2 Pi Aei + ^ei^ei 

»=1 i=l 
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defines the Poisson mapping between (T*V n , r , {■, ■}) and (gl(n), Ai) which is invariant with 
respect to the right S l O(r)-action (|2.12[) and the transformation (X,P) H» (— X, — P). In- 
deed, O is the composition of the following two Poisson mapping: 

G 1 (X,P) = $(X,P)+X, (X,P)£T*V n , r , 

e 2 (£, Y) = € + yy t , £ e so{n), y e M„ >r (E). 

The mapping 0! realizes T*V n r as a coadjoint orbit in the dual space of the semi-direct 
product so(n)(BpM nyr (M.) (e.g., see equation (29.11), page 225, [2B]; here p denotes the usual 
multiplication of matrices) and 02 is a Poisson mapping between (so(n) (B p M ntr (J&.))* and 
(so(n) © ac j Sym(n))* (see Lemma 7.1 of [44 ). 

We have also that the algebra of integrals (|5.10|) is the pull-back of (|6.4j) : 

$ = <d*F. 

However, the image Q(T*V n>r ) is the union of singular symplectic leaves in (gl(n) , Ai) . 
Namely, the generic symplectic leaf in (gl(n), Ai) has the dimension n 2 — n, while the 
dimension of generic leaf in Q(T*V n , r ) is 



r(r — 1) 



r 



21 = 2r(n - r) + ' - - ^' - -J < n 2 - n (6.5) 



(see Lemma [6.21 below). Nevertheless, it can be proved that the set of the functions 
is complete on a generic orbit laying in 0(T*V^ )r ) as well (see Lemma |6~51 below) . That is, 
among the integrals J- there is at least I polynomials p\ , . . . , pi independent on the symplectic 
leaves in 6(T*F„, r ). 

The rest of the proof follows the idea of [HJ [51]. Namely, since SO(r) acts on T*V n ,r 
freely and preserves the Poisson bracket {-,-}, the quotient space (T*V n ,r)/SO(r) carries 
natural induced Poisson bracket {•, •}'. Let a : T*V n , r — > (T*V n ,r)/SO(r) be the canonical 
projection. Then, by definition, 

{f,g}'(a(X,P)) = {F,G}(X,P), F = f o a, G = goo. (6.6) 

The Casimir functions ji, . . . ,j[ r /2] of the brackets {•, •}' can be obtained from the SO(r)- 
invariant functions Jk = tr(^> 2k ) via jk ° <J = Jk- 
The mapping induces Z2-Poisson covering 

: ((T*V n , r )/SO(r), {; •}') -> (e(T*F n>r ), A x ), d o a = 6. 

Hence the functions pioQ, . . .pioO are independent on a generic symplectic leaf and, together 
with the Casimir functions ji , . . . , j[ r /2] > form a complete commutative set of functions within 
(T*V n . r ) / SO(r). In other words, the functions 

Pioe,...,pjo6, Ji,...,J^j (6.7) 

form a complete commutative set of functions in the algebra of S'0(r)-invariant functions on 
T*V n ,r- Note that the independency of the functions Ji, . . . , J[ r /2] at (X, P) is equivalent to 
the regularity of Adgo( r )-orbit of ^(X, P). Therefore, (|6.1[) holds for the invariant manifolds 
where the rank of ^(X, P) is maximal. 

Then, according to Theorem 1 in [T3], the functions (|6.7[) together with \E'*(C 00 (so(r)) 
form a complete non-commutative set of functions in T*V nir . Therefore the integrals 0*^+ 
ty* (C°° (so(r)) of the Neumann system form a complete non-commutative set. Moreover, 
the functions (|6.7I) commute with all the integrals and therefore their Hamiltonian flows 
generate generic leaves of the isotropic foliation given by 0*J r + v I'*(C 00 (so(r)). Hence, the 
dimension of the generic isotropic tori is given by (|6.1[) . □ 
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Remark 6.1 The Hamiltonian of the Neumann system with the normal metric has the 



form H neumo = H o 9, where 



H = ±(h,h)-(A,v), 



defines the completely integrable (by means of the integrals (|6.4p ) Hamiltonian flow 

h=[v,A], v=[h,v] (6.8) 

on (so(n)®adSym(n))* . This system belongs to the class of Clebsch-Perelomov-Bogoyavlenski 
rigid body systems ([Ml [8], see a l so Section^]). Note that the representation of the Neumann 
system on the sphere as a system on a adjoint orbit is given by Ratiu |42) . 

Lemma 6.2 The dimension of a generic symplectic leaf in (Q(T*V n . r ), Ai) is given by 
formula \6.5\ ). 



Proof. Without loss of generality, choose a point Xq = (Ei, . . . ,E r ) S V n>r and a generic 
point (X 0) P ) € rj o V n ,r- Denote 

h = $(X ,P ), v = X Q x£ 

Then h is a generic so(n)-matrix of the form 

h = f h ^ , hi G so(r), h 2 G M r>n - r (M.) and v = ^ 

^From the definition (j6.3[) of the tensors Ai, A2 we get £ + 7; G kcr Ai(h + v), £ G so(ra), 
77 G Sym(n) if and only if 

M=0, (6.9) 
[t,h] + [r),v]=0. (6.10) 

The first equation gives the condition for £ to belongs to the subalgebra so(r) © so(n — r) . 
Denote 



67' ^"W ?72 



where £1 G so(r), £2 G so(n — r), 771 G Sym(r), 772 G Sym(n — r), 773 G M rj „_ r (R). Since 



-m 
4 

from (|6.2|) and (|6.10j) we find that = 0,, £2, 771, 772 are arbitrary, and 773 is uniquely 

determined from the equation 



6 o\ / 



'2 



bl'K-hl 



-7/3 

4 



= 0. 



For generic hi, the solutions of = form a maximal commutative subalgebra of 

so(r). That is, the dimension of the space of the solution of (|6.9|) . (|6.10j) is 

dim kcr Ai(h + v) = rank so(r) + dim so(n — r) + dim Sym(r) + dim Sym(n — r). 

Finally, the dimension of a generic symplectic leaf in (Q(T*V n , r ), Ai) is equal to n 2 — 
dimkerAi(/i + 7j) = 21. □ 
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Lemma 6.3 If all eigenvalues of A are distinct, then the set of the functions J6'.^| ) is a 
complete commutative set on a generic symplectic leaf in (0(T*V^ if -), Ai) 

Proof. We keep the notation from the proof of Lemma 16.21 The proof presented here is a 
modification of that of Theorem 1.6 in [5]. According to Theorem 1.1 in [5], the set of the 
functions (|6.4[) is complete on the symplectic leaf containing the point x — h + v if and only 
if 

(Al) All the brackets A.\ lt \ 2 non-proportional to Ai have the maximal rank n 2 — n at x. 

(A2) The kernel of the bracket A^-i at x, restricted to the linear space ker Ai, has dimen- 
sion n, 

dim{(£ + 77 G kerAi(x) | Ai,_i(£ + r?,ker Ai(x))| x = 0}. (6.11) 

Here all the objects are assumed to be complexified. Since the conditions (Al) and (A2) 
are both generic, it is sufficient to find x\ G Q(T*V n ,r) for which (Al) holds and X2 G 
Q(T*V n ,r) which satisfies (A2). Then the set of x G Q(T*V n . r ) satisfying both conditions 
will be open and dense everywhere in the induced topology on &(T*V nir ). 

If Ai + A2 7^ 0, A2 7^ 0, then the bracket Aa 1 ,a 2 is isomorphic to the canonical Lie-Poisson 
bracket on gl(n, C). Thus the brackets A\ lt \ 2 (Ai + A 2 7^ 0, A2 ^ 0) have the maximal rank 
n 2 — n at x = h + v if and only if the complex line 

C = {h + v - XA I A G C} 

intersects the set of singular points of gl(n, C) only at x = h + v. This condition is obviously 
satisfied if all the eigenvalues of A are distinct. 

To prove (Al) we have to find the (complex) dimension of ker Ai._i. From the definition 
(UTS) we get 

Ai,-i(£ + 77, C + 0)\h+v = ~(h, [v, 0}) + (A, [£, 9} + [ V , C]> (6.12) 
and ^ + 77 G kerAi + v), £ G so(n, C), 77 G Sym(n, C) if and only if 

[rt,A]=0, (6.13) 
[f J A]-[» 7 ^] = 0. (6.14) 

The solutions of 13[) are all diagonal matrices. For the given diagonal matrix 77, the matrix 
£ is uniquely determined from (I6.14[) . Therefore, dimker Ai i _i(/i + v). 

It remains to check the condition (A2). Take h of the form 

h= ( k o 0)' 

hi being a generic clement of so(r). From the proof of Lemma 16.21 we have 

ker Ai(ft + v) = so(r, C) hl ® Sym(r, C) © gl{n - r, C) , (6.15) 

where so(r, C)/ 11 = {£1 G so(r,C)\ = 0} is a Cartan subalgebra of so(r, C). 

^From (|6.12p . (16. 15)) we conclude that £ + 77 belongs to ker Ai._i(x)| kor Ai(z) if an d only 

if 

[Ai,ti]-[hi,rn]=0, (6.16) 

L4 llf7l ] Gso(r,C) x , (6.17) 

6 = 0, (6.18) 

[^2,772] = 0, (6.19) 
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where 771 and 772 are denned in the proof of Lemma 16.21 and 

A 1 = diag(a x , . . . , a r ), A 2 = diag(a r+1 , . . . , a„). 

Equations (I6.16|) . (16.171) form a closed system within g/(r, C), and from the proof of 
Theorem 1.6 [5], the dimension of solution of (|6.16[) . (|6.17l) is equal to r. On the other side, 
the solution of (|6.19[) are all diagonal matrices in gl(n — r, C). Hence (|6.1ip holds. □ 

Singular matrices A and the Einstein metrics. Now suppose that not all the eigen- 
values of A are distinct: 

a% = ■ ■ ■ = dfcj , aj. 1+1 = • • • = ak 1 +k 2 > • • • > a n+i-k r = ■ ■ ■ = a n , k% + k 2 + ■ ■ ■ + k r . 

Then we have a non-trivial isotropy algebra 

so(n) A = {£ G so(n) | [£, A] = 0} = so(ki) © so(k 2 ) © • • • © so(k r ). 

Let Q be the set of linear functions on so(u)a- The set T-\-Q is a complete non-commutative 
set of function on (gl(n), Ai) (Theorem 1.5 [§])■ By modifying the proof of Lemma 16.21 and 
Bolsinov's Theorem 1.5 [5], one can prove that the set of functions T + Q is complete on 
®{T*V n ^ r ) as well, implying non-commutative integrability of the Neumann systems (|5.3p . 
The complete verification is out the scope of this paper. 

Let SO{ti)a — SO{k{) x SO(k 2 ) x • • • x SO(k r ) C SO(n) be the adjoint isotropy group 
of A. The momentum mapping of the left 5'0(n) J 4-action is given by 

$>A = pr so(n)A $ = Wso(n) A (PX T - XP T ) 

and Q*Q are exactly Noether integrals arising from the S'0(n) J 4-symmetry of the Neumann 
flows. 

In particular, when A = 0, we get the integrals of the geodesic flow of the metric dsf. in 
the form $*(C°°(so(n)) + ** (C°° (so(r)) + where now 

S = {tr(A(PX T - XP T ) + XX T ) k I k = 1, . . . , n, A G K}. (6.20) 

We shall mention the following important corollary of the above construction. 
Let B : so(r) — > so(r) be positive definite and hs = |(£,B£). Suppose that the Euler 
equations 

/ = {/,Mr<,(r) e=K,vft B (0] = [e,Be], ee^w (6.21) 

are completely integrable with a complete commutative set of functions B. 

Theorem 6.4 (i) 03 + J? is a complete commutative set of SO (n) -invariant functions on 
T*V n , r; where % is given by WJU\) and 03 = **(£). 

(ii) XTie geodesic flow of the SO(n) -invariant metric ds^ on V n . r defined by the Hamiltonian 
function 

H n = - tr(P T P) - - tr ((X T P - P T X)M(X T P - P T X)) 

is completely integrable in the non-commutative sense. The complete set of first inte- 
grals is $*(C°°(so(n)) +5 + 23, and the generic trajectories of the system are quasi- 
periodic over the isotropic tori of dimension 

ddim(3 + 03) = dim V n>r - I. 

Here 21 is the dimension of a generic adjoint orbit in <$>(T*Vn,r) (see Lemma \2.5\) . 
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For example, if the matrix B = diag(i?i, . . . ,B r ) has distinct eigenvalues, we can take 
B(£) = B£ + £B (the Manakov operator on so(r)) and the commutative set 

23 = { tr(* + /iS) fe , n e E, k = 1, 2, . . . , r }. (6.22) 

Now, by using the chain method for the construction of complete commuting sets of func- 
tions on Lie algebras developed by Mikityuk [34] , one can easily prove that the Hamiltonian 
function (|3.14j) defines completely integrable system on so(r). Therefore we get 

Corollary 6.5 The SO(n) -invariant geodesic flows determined by Hamiltonian functions 
iS.13\) are completely integrable. In particular, the geodesic flow of the Einstein metrics 
constructed in 01/ are completely integrable. 

Note that the set of functions given in Theorem 16.41 differs from those described in 
Theorem 2.6. Another proof of integrability of geodesic flows of the Einstein metrics, based 
on using singular Manakov flows, is recently obtained in |17j . 

Commutative integrability of the Neumann flows. We turn back to the Neumann 
flows. According to Theorem 16. 1[ systems (|5.3p are integrable in the noncommutative sense 
by means of the integrals fy*(C co (so(r)) + where J is given by (|5.10p . However, the 
Neumann flows (15.3[) are integrable in the commutative (Liouville) sense as well: instead 
of (C°° (so(r))) one should take, for example, the commutative set (|6.22j) . Moreover, it 
follows: 

Corollary 6.6 Let all the eigenvalues of A be distinct. Suppose that the Euler equations 
16.21]) are completely integrable with a complete commutative set of functions B. Then the 
Neumann system with the kinetic energy given by the SOin) -invariant metric ds^ 

H neum , M = \ tr(P T P) - i tr {{X T P - P T X)M(X T P - P T X)) + ^ tr(X T AX) 

is completely integrable. The complete commutative set of first integrals is 5 + 23, where $ 
is given by HUty and 23 = **(#). 

In particular, the Neumann systems with the kinetic energy determined by SO(n) -invariant 
Einstein metrics constructed in 14-1 are completely integrable. 

7 Reduction to Grassmannians 

By definition, the points of the oriented Grassmannian variety G n ^ r are r-dimensional ori- 
ented planes passing through the origin in the Euclidean space M. n . The usual action of the 
group SO(n) on R n yields a transitive action on the set of all r-dimensional planes, i.e., on 
G n , r - The isotropy group of the r-plane spanned by the vectors E\, . . . , E r (relative to the 
base (|2.ip ) has the form 

{ SO r) SO(n-r))^ S0 ^* S °^-^ 

It follows that G n , r = SO{n)/{SO{r) x SO{n - r)). 

The oriented Grassmannian can also be seen as a quotient space of the Stiefel manifold 
by the right 5'0(r)-action described in Section 1. The quotient mapping V n _ r — > G n r is 

X = (ei, . . . , e r ) i — > e\ A • • • A e r . 
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The symplectic leaves in (T*V n , r )/SO(r), {■, •}') with {■, ■}' given by |0|} are the Marsden- 
Weinstein symplectic reduced spaces of T*V n . r . In particular, the reduced space that cor- 
responds to zero value of the momentum mapping 

*-\0)/SO(r), 

is symplectomorphic to the cotangent bundle T*G n ^ r equipped with a canonical symplectic 
structure. Note that (X, P) belongs to * _1 (0) if and only if X T P = 0. 

The last condition also implies that, although the Neumann systems (|5.4p and (|5.5|) are 
different on the whole T*V n ^ r , they coincide on ^^(O), hence their reductions onto the 
cotangent bundle T*G n , r are the same. 

The reduced flow can be written in the alternative Euler-Lagrange form. Namely, we 
have 

d r 

— (e! A ■ ■ ■ A e r ) = ei A ■ ■ ■ A e- t A ■ ■ ■ A e r 

and, in view of the matrix equation (|5.5[) . the reduced system is 

d 2 r 

(ei A • • • A e r ) = - ^ ei A • • • A Ae, A • • • A e r 

i—l 

+2 V" ex A • • • A e.i A ■ • • A e 3 A • • • A e r + A (ei A ■ ■ • A e r ) , (7.1) 

l__*<J__ r 

r 

A = tr A = tr{X T AX - X T X) = ^ ((e,, Ae t ) - (e 4 , e 4 ), ) . 

i=l 

where 

(ei,ej) = Sij, (e t ,e 3 )=0, i,j = l,...,r. (7.2) 

Note that (|7.ip is 50(r)-invariant. We refer to ()7. 1[) as a Neumann system on the 
oriented Grassmannian variety G n ,r- 

Theorem 7.1 Suppose that all the eigenvalues of A are distinct. Then the Neumann system 
on T*G rhr is completely integrable in the Liouville sense by means of the integrals induced 
from the SO(r) -invariant functions \5.10}) . 

Proof. Since T*G n , r = f- 1 (0)/SO(r) is a singular symplectic leaf in T*V n , r /SO(r), the 
statement of the theorem does not follow directly from Theorem 16. II 

We keep the notation from the proofs of Theorem 16.11 Lemma 16.21 and Lemma 16.31 Let 
h = ty(X ,P ), v = X Q X£, where (X ,Pq) S * _1 (0) n T^ Q V n . r is in a generic position. 
Then hi = 0. From the proof of Lemma l6.2l we get 

dimker Ai(/i + v) = dmi(gl(r) ® gl{n — r)). 

and the mappings and map V I' _1 (0) and 5 ,_1 (0)/S'O(?*) to the single symplectic leaf 
in (gl(n), Ai). Now, by modifying the proof of Lemma 16.21 one can prove conditions (Al) 
and (A2) for x = h + v with hi = as well. Therefore, among functions (|6.4p one can find 
exactly 

r(n — r) = — (dimgZ(n) — dim^(r) — dimgl(n — r)) = dimG Iljr 

independent functions Pi,. ■■ ,p r (n-r)- Thus the functions p\ o 9, . . . ,p r („_ r ) ° 8 provide a 
complete commutative set of functions on ^/ _1 (0)/S'O(r). □ 
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The special case r — 1. The Stiefel variety V n , n -x is diffeomorphic to SO(n) = V n y. to 
ei, . . . , e n _i one can associate the unique unit vector e n such that ei, . . . , e„ is the orthonor- 
mal base with the same orientation as E\,..., E n : 

X = (ex, ■ • • , e„_x) i — S> = (ex, . . . , e n _x, e„) £ S"0(n). (7.3) 

Similarly, the oriented Grassmannian variety G„. n _i is diffeomorphic to G n _\ — S n via 
the mapping 

d A • • • A e n _i i — > e„. (7.4) 

It is natural to expect that in this case the Neumann system (|7.ip gives rise to the 
classical Neumann system on the sphere Indeed, due to conditions (17.21) . for r — 1 the 

second term in the right-hand side of equations (|T. 1 [) vanishes. Then, under identification 
(|7.4p . these equations gives rise to 

e„ = (A- tiAI)e n + Ae„ = Ae n - ((e n ,Ae n ) + (e„,e„)) e„ , (7.5) 

which describes the motion on the sphere 5 n_1 = {(e„,e„) = 1} with the potential 

The 4-th degree potential. By using Cartan models of symmetric spaces, a class of 
new integrable potential systems on such spaces was obtained by Saksida [45j . As noticed 
in in the case of the sphere 5 1 ™ -1 , such a system is a generalization of the Neumann 
system in presence of a potential of degree 4. The latter system is separable in the spherical 
elliptic coordinates and was found previously in |50) . 

In addition, it can be proved that the construction of |45j on the Grassmannian varieties 
gives rise to the Hamiltonian 

H(X, P) = 2 tT ( pTp ) + tr(X T A 2 X) - tr(X T AXX T AX), (7.6) 

which for r = 1, takes the well known form H = \(p,p) + J2i=i a 1 e ? — a i e T) 2 [501146]. 

For r > 1, the Hamiltonian flow on T*V n _ r determined by the Hamiltonian function (|7.6[) is 
integrable after its restriction to the invariant manifold \E' _1 (0) C T*V n . r and the reduction 
to T*G ntr . This system will be discussed elsewhere. 



8 The Dual Lax Pair and Geometric Interpretation of 
Integrals 

As mentioned in Section 5, like the classical Neumann system on the sphere S n ~ 1 , the 
Neumann systems on V(r, n) also admit dual Lax representations. 



Theorem 8.1 Up to the action of a discrete group generated by reflections $4-(ty , the 
Neumann flows |5.^| ) and \5. 5\) are equivalent to the following 2r x 2r matrix Lax pair with 
a rational spectral parameter A 

"TL^neum^V = l^netim (^) ' -^neura (^)] > (^-l) 

r * ( -x T (\i n -A)-ip -X T (XI n — A)~ X X\ 

A-neun. W - I ^ + p^(AI„ - A)" X P P T (AI„ - A) -1 ^ J ' ^ J 
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where for system \5.J$ , respectively H5.5]) , one should put 



( X T P I \ / I \ 

-4ne«™( A ) = ( A _ AI? _ pT x J . respectively, .4£ ettm (A) = ( A - AI r ) ' ( ' 8 ' 3 ' ) 

w/iere A = X T AX - P T P. 

The statement is checked straightforwardly by using constraints (12. 4p and the matrix 
identities 

A{\I n - A)^ 1 = (AI„ - A)- 1 ^ = A(AI„ - A)- 1 - I„ . 



The dual Lax pair ()8.3p for the generalized Neumann system (|5.5[) was first given in 
unpublished manuscript |28j . For r = 1 it gives the known 2x2 Lax pair for the Neumann 
system indicated in several publications (see, e.g., [47] and references therein). 

Remark 8.1 For r the considered Neumann flows become geodesic flows of bi- invariant 
metrics on SO(n). Then both Lax representations (n x n and 2n x 2n) give the integrals 
for the geodesic flow of a bi-invariant metric, but also they give the integrals for the right- 
invariant Manakov geodesic flows on SO(n) described by equations Q3.5p . (I4.2[) . 

The spectral curve and integrals. Let 

a(A) = (A - oi) • •• (A - o n ). 

The spectral curve of jC-^ eum (X) can be written in form 

\a(\)C* neum (A) - win | = w 2r +w 2r - 2 a(\)2 2 {\) + - ■ • + W 2 a 2 ''- 3 (A)I 2r _ 2 (A) + a 2r - 1 Z 2r (A) = 0, 

where I21 (A) are invariant polynomials in the components of the wedge products 
ej 1 A • • • A ej. and 

ei A • • • A e r , 

e\ A • ■ • A e r A pi, i = 1, . . . , r, 

(8.4) 

ei A • • • A e r A pi A • • • A p r . 

Note that, due to the symplectic block structure of £* eum (A), the coefficients at odd powers 
of w in the spectral curve are zero. 
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In the case 2r < n the polynomials can be written in form 
Z 2 (A) 



where Ik = {ii, ■ ■ ■ ifc} C {1, . . . , n} is the multi-index with distinct indices 1 < i\ < • • • < 
ik < n and l^l/* is the k x k diagonal minor of the momentum matrix $ corresponding to 
the multi-index Ik- Note that, in view of definition of $, 

\®\iZ = (ei A • • • A e r A pi A • • • A p r )j 2r - 

In the case 2r > n the polynomials l2/(A) have the same form with the only difference: 
the terms with the wedge products of €i,pj of order > n are absent. 

It follows that in both cases X 2 i{\) are polynomials in A of degree n — I and that the 
leading coefficients of 2a (A), . . . ,Z2 r (A) produce trivial constants on V(r,n). Hence, as a 
simple counting shows, the number of nontrivial integrals on T*V(r, n) provided by the Lax 
matrix £* emn (A) in (|8.3p equals 

N = (n - 1) + (n - 2) H h (n - r) = r(n - r) + r(r - l)/2, 

which coincides with the dimension of the Stiefel variety. 

Note that, although the Lax matrix £* eum (A) is not invariant under the right SO(r)- 
action, the spectral curve and therefore all the integrals l2i(A) are S'0(r)-invariant. 

Since the number TV is bigger than half of dimension (|6.5I) of a generic symplectic leaf 
within (T*V n . r )/ SO(r), some of the integrals are dependent. Like the "big" Lax matrix 
£neum(A) in (|5.7|) . the dual Lax matrix £* eum (A) does not produce explicitly the momenta 
integrals ^ij. 

Geometric interpretation of the integrals l2r(A). The components of forms (|8.4p that 
appear in the last invariant polynomial l2r(A) have a transparent geometric interpretation: 
they are Phicker coordinates of the 2r-dimensional linear subspace (2r-plane) 

£ = S(X, P) C R n+r (x u . ..,x n , yi ,..., y r ) 
24 



h Ji=\ji,—,n} 

+ E (A _ a . ) a( . A ( ) A _ a . ) E (e A A...A ejl / 
I l+1 {A a ^> [A ai ^> Jt,l<i<r 

= E n T^Tx \ < ei A ' ' ' A e ^ 

+ Ea 7 E( ei A-'-AerApi)?, 

J r + 1 * — 1 

+ Et7 7^T7 7 E (eiA---Ae r A 

+ v ^ 1$ Ii 2 " 

+ (A-a n )...(A-a 42 J 1 



spanned by the columns of the 2r x [n + r) matrix 



V = V(X,P) = 



/ei ■•• e r px ■■■ p r \ 

••• 

: 1, 

\0 •■• J 



(8.6) 



I r being the identity r x r matrix. Indeed (see, e.g., [25] ). for any k,m(k < m), the Pliicker 
coordinates of a fc-plane 7r in W a (xi, . . . , x m ) spanned by independent vectors v±, . . . , v k £ 
R m are the coefficients Gi of the polynomial 



vi A • • • A v k = ^2 Gi dx ix A ••• A dx ik , 



where / = . . . ik} C {1, . . . , n} is the multi-index with 1 < i\ < ■ • • < < n. 

Then the Pliicker coordinates of E are given by all 2r x 2r minors of V. In particular, 
the 2r x 2r minors that completely contan I r give the Pliicker coordinates of the r-planc 
span(ei, . . . , e r ) C W l . 

Now consider the following family of confocal cones in R ra+r (a;i, . . . , x n , y\ . . . , y r ) 

Q(X) = \ Y ^ + --- + T ^ + yf + --- + y 2 r=o), AgR. (8.7) 

I.A- 01 A — a n J 

The following theorem is a first variant of a generalization of the remarkable Chasles the- 
orem describing a geometric relation between the geodesic flow on an ellipsoid and common 
tangent lines of confocal quadrics f [l5l l35l |30] V 

Theorem 8.2 Let the 2r-plane E(i) C R n+r be associated to a generic solution (X(t),P(t)) 
of the Neumann systems |5.^[ j and i5.5\) on V n , r as described above. Then E(i) is tangent 
simultaneously to n — r fixed confocal cones Q{c\) 1 . . . , Q(c n _ r ), where ci, . . . , c„_ r are the 
roots of the invariant polynomial l2r(A). 

One can show that for real solutions (X(t), P(t)) all these cones are real. 

In the particular case r = 1, one can also consider the section of E and of the family Q(X) 
by the subspace {y± = 1} = M™, which give respectively an afhne line l(t) = p(t) + span{e(£)} 
and the family of confocal quadrics 



^ a\ — a a n — A 

Then, due the above theorem, l(t) is tangent to n — 1 fixed quadrics Q{c\), . . . , Q(c n _ r V, 
and we recover the following variant of the Chasles theorem given for the Neumann systernj 

Proposition 8.3 (Moser, [55] ) 
(i) Let (e(t),p(t)) be a solution of the system on r p*S n ~ 1 with the Hamiltonian 

H = ±a i F i , ^^jirf, 
^ — ' ^ — ' <n — a, 

cti being arbitrary constants. Then the associated line l(t) is tangent to n — 1 fixed 
confocal quadrics of the family Q(X). 



1 lt can be formulated in two different ways (see Theorem 12 in |20| given for the Clebsch-Perelomov 
systems and Theorem 4.10 in [37]). 
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(ii) If (e(t),p(t)) is a solution of the system on x*S n ~ 1 with H — X)™=i-^V a * restricted 
to W = 0, the corresponding line l(t) is tangent to the ellipsoid Q(0), on which the 
contact point I D Q(0) traces a geodesic. 

The proof of Theorem 18.21 is based upon the following property described in [5D] . 

Proposition 8.4 Let Gi, I — {ii, . . . ife} be the Pliicker coordinates of a k-plane passing 
through the origin in R m . The set of all such k-planes that are tangent to a nondegen- 
erate cone {(x,Bx) = 0}, B = diag(6i, . . . b m ) is the intersection of the (non- oriented) 
Grassmannian G(m, k) C A fc R m with the quadric 



5>ijg? = qJc 



A fe R"\ \B\j = b n ---b lk . (8.8) 

Proof of Theorem 18.21 Let us now set m + r and consider the family of cones (|8.7|) . Let, 
as above, £ be a 2r-plane in R™ +r spanned by the columns of V and associated to a point 
X,P on T*V n>r , and let 

( 1 1 
B = diag , . . . , , 1, . . . , 1 

\ A — a\ A — a n 

Then, in view of Proposition 18.41 and the structure of V, for a fixed A = A* the set of the 
2r-planes that are tangent to Q(X*) is described by the following quadratic equation in 
terms of the Pliicker coordinates of £ 

Ir (A* - a h ) - (A- a ir ) + g (A* - a*) • • • (A* — a w ) g^ 1 

+ £ (A*- an )..V- a , ~) E (eiA-Ae r A A Ap i )?^+- 

/ r . +2 v iy v r+2; l<i<j<n 

+ ^ (A-- atl ).. 1 (A*- ai2 ) (eiA --- Ae - ApiA --- A ^ =0 - 

Due to (I8.5p . this coincides with the equation X2 r (A*) = up to multiplication by a(A*). 

Since X2 r (X) is also an invariant polynomial of degree n — r in A, for a fixed plane £ 
there exist precisely n — r fixed cones of the family Q(\) tangent to £. This establishes the 
theorem. □ 

Restriction to R". By analogy with the case r = 1, one can consider the restriction of 
family (18.7[) to the linear subspace {yi = ■ ■ ■ = y r = 1}: 

QrW = i^K + ■■■ + — =r}= i" 1 (Q(A) n { yi = 1, . . . ,y r = 1}) , (8.9) 
t ai — A a n - A J 

where i : 1" — >• R ra+r is the natural inclusion 

i(xi , . . . , Xji) — (^i , . . . , x n: 1, . . . , 1). 

This gives a family of confocal quadrics in R™. 

Further, the section of £ by the subspace {y\ = • • - = y r = 1} defines an affine r-plane 
£(t) = i _1 (£ n {y\ = ■ ■ ■ = y r = 1}) in R n (xi, . . . , x n ), which is spanned by the orthogonal 
vectors ei, • • • , e r and passes through the point p% + • • • + p r - As a result, to a generic 
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solution (X(t), P(t)) of the Neumann system on V n ^ r one can uniquely associate the moving 
r-plane 

E(t)=p 1 (t) + ---+p r (t)+span{e 1 (t),--- ,e r (t)}. 

In contrast to the case r = 1, due to dimensional reasons, for r > 1 the r-plane £(i) is not 
necessarily tangent to the quadrics Q r {c\) 1 • • • , Q r {c n -r)- More precisely, since 

di (T( Xu ...,x n )QW) = Ti(x u ...,x n )QW n{y% = l,...,y r = 1}. 

the tangency of £(i) and Q(ci), for a fixed t, either implies the tangency of the corresponding 
affine r-plane £(i) and the quadric Q r {ci), or £(i) does not intersect Q r (ci). As a result, 
one cannot formulate a natural generalization of the Chasles theorem in R ra that involves 
this r-plane. 

Another feature of the case r > 1 is that, although the first integrals given by the poly- 
nomial l2r(X) are invariant with respect to the right S'0(r)-action on (X, P), the 2r-plane £ 
and r-plane £ do not have this property. Thus, a generic polynomial l2 r {\) corresponds to a 
whole family of 2r-planes (r-planes, respectively) that are tangent to the same set of confocal 
cones and is obtained as the orbit of £ (£, respectively) under the right S'0(r)-action. 

Then, it natural to replace £ by the moving cylinder A(i), the union of 2r-planes 
Y,(X(t)B, P(t)B) spanned by the columns of the 2r x (n + r) matrices 

V(X(t)B,P(t)B), BeSO(r), 

where V(X,P) is given by (|8.6p . The cylinder A(t) is 50(r)-invariant and, due to the 
construction, is tangent simultaneously ton-r fixed confocal cones Q(ci), . . . , Q(c„_ r ). 

Next, the section of A(t) by the subspace {y\ = ■ ■ ■ = y r = 1} defines the moving 
(2r — l)-dimcnsional cylinder 

A (*) = jX^'i^)! 56 5 °( r ) j + span{ei(t),-- - ,e r (t)}, 

which is now an appropriate object for the second generalization of the Chasles theorem: 

Theorem 8.5 Let the (2r — 1) -dimensional cylinder A(t) C K n be associated to a generic 
solution (X(t), P(t)) of the Neumann systems J5.^[ ) or \5.5\) on V n ^ r as described above. 
Then A.(t) is tangent simultaneously to n — r fixed confocal quadrics Q r (c±), . . . , Q r (c n - r ) 
of the confocal family J#. 7| ), where c\, . . . ,c„- r are the roots of the invariant polynomial 
l 2r {\). 

Proof. First, note that the plane E(X(t)B, P{t)B) can be obtained from E(X(t), P(t)) by 
rotating it in the coordinates yi,...,y r by the matrix B . That is, the cylinder A(t) can 
be regarded as the orbit of £(X(i), P(t)) with respect to the 5'0(r)-action in the coordinates 
yi, . . . ,y r . This property is related to the S'0(r)-symmetry of the cones Q(X) in y\, . . . , y r . 
Indeed, the 2r x (n + r) matrices 

ei • • ■ e r pi ■■■ p r \ 
••■ 

; ; b- 1 

■•• / 

define the same 2r-plane T,(X(t)B, P(t)B), whereas the second matrix is obtained from 
V(X, P) by left multiplication by the block matrix diag(I„_ r , B~ x ). 



V(XB, PB) and V(XB, PB) 
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Now let l(t) be the line along which the plane H(X(t),P(t)) is tangent to the cone Q(cj), 

I = l(t) = spanM*) = (v u v 2 , . . . , v n+r ) T } C £(X(i),P(i)) C R" +r 

and I sit) = span{uB(i)} be the tangency line of the rotated plane T,(X(t)B, P(t)B) and 
Q(ci). Due to the above observation, vs(t) = diag(I„, B~ 1 )v(t). 

One can always find B* £ SO(r) depending on i and t, such that the last r coordinates 
of VB*(t) are equal. Then the r-plane 

E B . (t) = r 1 (fl(X(t)B*,P(i)B*) n {yi = • • • = y, = 1}) = £ BJ>(t) + span{ ei , • • • , e r } 

is tangent to the quadric Q r {ci) at the point i _1 (£b*(£) PI {t/i = ■ ■ • = y r = 1})- Since for 
any i and i, £s*(t) is a subspace of the cylinder A(t), we arrive at the statement of the 
theorem. □ 



Chasles Theorem for Manakov Flows. Similar statement holds for the geodesic flows 
of submersion metrics defined by Manakov operators (see also Theorem 12 in |20j). 

Suppose 2r < n. The dual Lax pair for the flow with Manakov operator (|4.4I) given in 
Theorem 14.41 gives the set of commuting integrals J72z(A), the coefficients with term w 2r ~ 21 
in the expression \a(X)C^ nan (X) — iwl n |. In particular, 

^ (A - OiJ • • • (A - a l2 J 

is a polynomial of degree n — 2r in A. 

Consider the following family of confocal cones in M. n : 



Qo(A) = 




A e M. 



Repeating the arguments of Theorem 18.21 and using the fact that the Manakov geodesic 
flows with different choices of the matrix B in (|4.2j) are quasi-periodic motions over the 
same isotropic toric foliation of T*V n ^ r , we get: 

Theorem 8.6 Let the 2r-plane 

S(t) = span{ei(i), . . . , e r (t), Pl (t), . . . ,p r (t)} C E" 

be associated to a generic solution (X(t), P(t)) of the geodesic flow \3. 5\) given by the Man- 
akov operator [4-<fy - Then £(i) is tangent simultaneously to n — 2r fixed confocal cones 
Qo(ci), ■ • ■ , Qo(cn-2r), where a, . . . , c„_2r are the roots of the invariant polynomial j72r(A). 

Remark 8.2 Since all objects in theorems 18.51 and 18.61 are right 50(r)-invariant they are 
also valid for the Neumann system on the oriented Grassmannian variety G n ,r as well as for 
the geodesic flows on G n , r obtained by submersion from the Manakov flows. 



The case of the Poisson sphere. As an illustrative example, consider the case r = 1 
and Manakov operator (|4.2I) with B = — A^ 1 . Then the submersion metric on 5' Tl_1 takes 



ds2 = jA^J2 a ^- ( 8 - 10 ) 



the following form (see [14] 
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In the elliptic coordinates the metric (18.101) is of the Stackel type, and its geodesic flow 

e = (A^e, e)A- 1 P - (A~ 1 e,p)A~ 1 e, 
p = (A~ 1 e,p)A~ 1 p - (A~ 1 p,p)A~ 1 e 

is completely integrable [14 . For n — 3, the metric (|8.10p is proportional to the metric on 
the Poisson sphere S 2 , i.e., to the metric obtained after SO(2) reduction of the free rigid 
body motion around a fixed point with the inertia tensor / = A . 
From Theorem l8.6l we get 

Corollary 8.7 Let e(t) be a geodesic line of the Poisson sphere metric \8.10)) . Then the 
moving 2-plane 

£(t) = span{e(£), p(t)} = span{e(i), Ae(t)} 
is tangent to n — 2 fixed confocal quadrics of the family Qq(\). 

9 Appendix 1. Magnetic Neumann Systems 

The Neumann system on S 1 ™ -1 in presence of the Yang Mills fields. We now go 

back to the reduction of the Neumann system to Grassmannians described in Section [7] and 
consider in detail the case r — 1. 

Apart from the complete reduction of the Neumann system onto '&~ 1 (0)/SO(n — 1) = 
T*^' 1-1 , it is also convenient to describe partially reduced flows for non-zero values of the 
momentum 'J, that is, the flows on 

(T*V n . n ^)/SO(n - 1) S so{n) x S n ~\ (9.1) 

as well as the reduced flows on the symplectic leaves 

*-\ri)/SO{n - 1), S V-^OJ/SOin - 1) C {T*V n , n ^)/SO(n - 1), 

where O v is the adjoint orbit of n. It is known that for n ^ 0, the quotients i ^~ 1 (O ri ) / SO(n— 
1) are diffeomorphic to 0,,-bundles over (co)tangent bundle of the sphere S n ~ 1 and that the 
reduced systems are natural mechanical systems with the influence of the Yang-Mills fields 
(see, e.g., [35], Chapter III, which provides a detailed geometrical analysis of such systems). 

We shall derive the reduced equations and Yang-Mills fields directly from the equations 
of motion. Namely, for a matrix £ 6 so(n) and an unit vector e n G S*" -1 , define projections 

P f e „ (0 = £ e « ® e„ 4- e„ <g> e„£, pr^x (£) = £ - pr e „ (f) 

with respect to the orthogonal decomposition 

so(n) = {e n A R™} © {e n A M"} 1 . (9.2) 

Note that {e n A R™} 1 - = so(n — 1) and {e„ A R™} can be naturally identified with the 
tangent space T eji S" l_1 . Therefore the reduced phase space T*V n , n -i/ SO(n — 1) is also 
represented as a so(n — l)-bundle over T5' n_1 that we shall denote by so(n — 1) x s TS n ~ 1 . 
There is the natural inclusion 

so(n - 1) x s TS 71 - 1 C so{n) x TS 11 - 1 : (£,e„,e„) G so(n - 1) x s ^ pr e „(£) = 0. 

Proposition 9.1 (i) The reduced Neumann system on the quotient variety \9. 1)) has the 

form 

Q=[A, e„®e„], e n = Qe n , Q G so(n). (9.3) 
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(ii) The second derivative of the vector e n is 

e n = Ae n - ((e n ,Ae n ) + (e„, e„)) e„ + F n e„, (9.4) 

where Fq = pr e?i _L (ft) = ft — &(e n , e n ) and $(e„, e„) = e n A e n is i/ie standard SO(n)- 
momentum mapping on TS 11 ^ 1 . 

If we restrict the flow to the invariant submanifold 4 ,_1 (0^)/S'0(r), then (i*b,e n ,e„) 
ranges over the subbundle O v x s TS n ~ 1 C so(n — 1) X s TS 71 " 1 obtained by replacing fibers 
so(n — 1) by the adjoint orbits C so(n — I). In particular, for n = we recover (|7.5[) . 
For ?y ^ 0, the additional term Fq e n can be interpreted as the influence of the Yang-Mills 
field with the internal symmetry group SO(n — 1) and charge type O v ( 26J). In the special 
case n = 3, r = 2, after identification of the Lie algebras (so(3), [•, •]) and (R 3 , x), equation 
takes the form 

e 3 = ^e 3 - ((e 3 , Ae 3 ) + (e 3 ,e 3 )) e 3 + ee 3 x e 3 , 

which describes the motion of the particle with the charge e on iS 2 in the magnetic monopolc 
field. Here 77 = eE\ A E2 and Fq e 3 = e e 3 x e 3 represents the Lorentz force of the magnetic 
monopole. 

Since we already proved the completeness of the commuting integrals (|5 . 10[) on a generic 
symplectic leaf within (|9.I[) , we arrive at the following statement. 

Theorem 9.2 The Neumann system perturbed by the Yang-Mills field, i.e, the restriction 
of i f 9. 3\) to O n x s TS n ~ l , is completely integrable for a generic value rj € so(n — I). 

Proof of Proposition [Ol Let us identify V n ,n-\ and SO in) via (|7.3p and consider SO(n) 
as the configuration space of the rigid body moving around a fixed point. Then the vectors 
ei, . . . , e„ are fixed in the body, the matrix X = (ei, . . . , e n ) maps the fixed frame to the 
frame attached to the body and 

M = VX T - XV T 

plays the role of the angular momentum of the body in the space frame. Here we denoted 
the n x n momentum matrix by V . 

From the identity e n C§) e n = I„ — ~Yld = i e, ® ei we see that, after identification (|7.3I) . the 
Neumann system with the normal metric on V"„. n _i corresponds to the motion of the rigid 
body with the Hamiltonian 

H C p = \{M,M) + i (to A -(t^, A^)) ■ (9.5) 

This is a special, symmetric case of the Clebsch-Perelomov rigid body problem: the 
inertia operator of the body is the identity on the Lie algebra so(n) (see [H]). Therefore the 
angular velocity in the space frame Q = X ■ X~ 1 and the angular velocity in the body frame 
lu = X~ x X are equal to the angular momentum in the space frame M and to the angular 
momentum in the body frame m = X~ X MX respectively. 

The Hamiltonian (|9.5j) is invariant with respect to rotations in R" _1 = span(ei, . . . , e n _i), 
i.e., it is right SO(n — I)-invariant. In the right trivialization, the motion of the body is 
described by the Euler-Poincare equations 

tl=[A 1 e n ®e n ] (9.6) 

together with the Poisson equations 

ej = Qej, i=l,...,n. (9.7) 
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Whence it is clear that the reduced flow is given by (I9.6[) and the last Poisson equation 

e n = fie n . (9.8) 

To prove the second assertion of Proposition 19. 1[ we rewrite equations (|9.3|) with respect 
to the orthogonal decomposition (|9.2p . Then, from (j9.8[) , we have e n = pr 6n (f2)e„ and 

pr en (Q) = e„ Ae„ = $(e„,e„). (9.9) 

In view of (|9.6[) and (|9.9j) . the time derivation of (|9.8j) reads 

e„ = Qe n + fte n = 0e„ + pr 6n (fi)e„ + pr 6n x (0)e„ 

= (Ae„ ® e n - e„ ® e„A)e„ + (e„ A e„)e„ + pr e ^ (Q)e„ 
= ^4e„ — ((e„, Ae„) 4- (e„, e„)) e n 4- Fq e n , 

which concludes the proof. □ 

Symmetric Clebsch Perelomov Bogoyavlenski rigid body systems. A similar re- 
duction can be made for the Neumann systems on V n>r and G n ,r by considering the motion 
of a symmetric rigid body with the Hamiltonian 

Hcpb = \ (M, M) + i (ei, Aa) , 

which corresponds to a special (symmetric) case of the Bogoyavlenski generalization of the 
Clebsch-Perelomov system (see [8]). Namely, in the space frame the motion is described by 
the Euler-Poincare equations 

f2=[ei®eH h e r <B> e r , A] (9.10) 

together with the Poisson equations (|9 . T[) . After substitutions 

h = Q and v — e\ ® e\ + • • • + e r ® e r 

they take the closed form (|6.8p . 

The system is right SO(r) x SO[n — r)-invariant with the momentum mapping 

* = *so(r) + *«,(n-r)i ^ so(r) = pr so(r) (A^Of), * so( „_ r) = pr so( „_ r) ( X^SIX) , 

where 

so(r) = span{Sj A Ej, 1 < i < j < r}, so(n — r) — span{.Ej A Ej, r + l<i<j< n}. 
The reductions of (|9T0]l . (|9~7| to *7 1 ( „_ r) (C7 ?so( „_ r) )/50(n - r) and to 

* -1 (0.?„ w x O nio(n _ r) )/50(r) x SO(n-r) 

lead to the Neumann systems on the Stiefel variety V n ^ r and, respectively, to the oriented 
Grassmannian variety G n ^ r under the influence of the Yang-Mills fields. In particular, if 
Vso(r) = (i]so(r) — Vso(n-r) = 0), we get the Neumann system with the normal metric on 
V rh r (respectively, the Neumann system on G„ jr ). 
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Magnetic Neumann flows on V n . 2 and G nt2 . The adjoint orbits in so(2) are points, so 
a symplectic reduced space ^ _1 ((j7 so (2)> 0))/S'O(2) x SO(n — 2) is diffeomorphic to T*G Ui2 . 
For rj so r 2 \ = ei?i A -E2 7^ 0, it represents the magnetic cotangent bundle T*G n , 2 - the 
canonical symplectic structure of T*G n , 2 is "twisted" by adding the magnetic form, which 
is exactly Kirillov-Konstant symplectic form on G n ^ 2 multiplied by e (for more details sec, 

e.g., mm)- 

As above, for us it is convenient to write the equations in the Euler-Lagrange form. Due 
to the presence of the magnetic form, the tangent bundle momentum mapping of the left 
S'0(n)-action is modified by adding the term ee\ A e 2 (see [T9"l 113j ) 

$ £ = ei A ei + e 2 A e 2 + e e\ A e 2 = [ei A e 2 , ei A e 2 + ei A e 2 ] + e ei A e 2 , (9.11) 

whereas the right-hand side of the Neumann system (|7.1[) is modified by adding the term 
e$o (see Q3]) 

^ (ei A e 2 ) = — Aei A e 2 - ei A Ae 2 

+ {{e x ,Aex) - (ei,ei) + (e 2 ,Ae 2 ) - (e 2 , e 2 )) ei A e 2 (9.12) 
+2ei A e 2 + e [ei A e 2 , ei A e 2 + e\ A e 2 ]. 

Here ei,e 2 ,ei,e2 satisfy the conditions (|7.2[) . (Clearly, a similar symplectic reduction with 
a magnetic term can also be applied on the Stiefel variety V n . n - 2 (see jH]).) 
To construct the magnetic Neumann flows on T*V n>2 , consider closed 2-form 



dei A de 2 = ^ de\ A de 



restricted to V n , 2 . Let 7r : T*V n>2 — > Ki,2 be the canonical projection and define the sym- 
plectic form 

ui e = u) + e n*L) mag , (9.13) 

where w is the canonical form on T*V n , 2 (see Section [2]). 

The following two propositions can be verified by straightforward computations. 

Proposition 9.3 The left SO(n)- action on (T*V n . 2 ,uJ(_) is Hamiltonian with the momen- 
tum mapping given by 

$ e = $ + e ei A e 2 = pi A e x + p 2 A e 2 + e e\ A e 2 . (9-14) 

Proposition 9.4 The Hamiltonian equations defined by the Hamiltonians of the Neumann 
systems with the Euclidean metric and the Normal metric with respect to the symplectic 
structure h9.13\) read 

ei =Pi, 

e 2 = p 2 , (9.15) 
Pi = -Ae\ + ((ei, Ae\) - - e{e\,p 2 )) e\ + ((ei,Ae 2 ) - (pi,p 2 )) e 2 + ep 2 , 

P2 = ~Ae 2 + ((e 1 ,Ae 2 ) - (px,p 2 )) e x + ((e 2 ,Ae 2 ) - (p 2 ,p 2 ) + e(e 2 ,pi)) e 2 - ep x 

and, respectively, 

ei = $ ei = pi - (ei,p 2 )e 2 , 

e 2 = $oe 2 =pi - (e 2 ,Px)e 1 , (9.16) 
Pi = $oPi - Ae x + ((e l5 Aei) - e(e 1; $ e2)) ei + {e u Ae 2 )e 2 + e $ e 2 , 
P2 = $oP2 - -4e 2 + (e 2 ,Ae 1 )e 1 + ((e 2 ,Ae 2 ) + e(e 2 ,$ ei)) e 2 - e<J> ei . 
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Equations (|9.15p . f|9 . are right SO(2)-invariant and have integral = (ei,p2) — 
(e2,Pi). The magnetic Neumann system f|9 . 1 2[) can be also seen as a reduction of the 
system (|9.15p . or (|9.16[) with respect to the right S'0(2)-action. 

Furthermore, for the systems (|9.12p . (I9.15|) . (|9.16l) the relation (jEB] still holds in the 
form 

= [ei (8) e 1 + e 2 ® e 2 , A], — (e x <g> e x + e 2 ® e 2 ) = [$ e ,ei ® e 1 + e 2 ® e 2 ], 

which implies the Lax representation (|5.7p , where instead of the momentum mapping $ one 
should use $ c in (|9.14l) . 

Theorem 9.5 27ie magnetic Neumann systems &9.12\) , \9.15\) and \9.16\) are completely 
integrable in the commutative sense with respect to the twisted symplectic structures described 
above. 

The proof is a simple modification of those of Theorems 16.11 and 17.11 

10 Appendix 2. Rank r Double, Coupled and Neumann 
Systems on Complex Stiefel Manifolds 

In this section we briefly consider several natural generalizations of the Neumann flows on 
Stiefel varieties. We present their equations of motion and Lax representations, however a 
complete verification of the integrability is out of the scope of this paper. 

Rank r double Neumann system. In [?7] Suris introduced the double Neumann system 
describing the motion of 2 points x, y S K™ which interact via the bilinear potential (x, Ay) /2 
under the constraint (x,y) = 1. 

We consider rank r double Neumann system defined by the Lagrangian function 

L(X, Y, X, Y) = tr(X T Y) - tr(X T AY), (10.1) 

where the n x r matrices X,Y £ M„ jT .(R) are subject to the constraints 

X T Y = I r . (10.2) 

The corresponding Euler-Lagrange equations with r x r matrix multipliers readH 

X = -AX + XA T , Y = -AY + YA, (10.3) 

where 

A = X T AY - X T Y. (10.4) 

The rank r double Neumann system is an extension of the Neumann system on V n , r 
with the Euclidean metric: if (X(t),Y(t)) is a solution of the system (|10.3|) with the initial 
conditions X — Y, X — Y, then (X(t),P(t)) = (X(t),X(t)) is a solution of ([53]). The Lax 
representation (|8.ip extends as follows. 

Theorem 10.1 The equations 111 0. 3\) imply the following 2r x 2r matrix Lax pair with the 
parameter A 

j t C(X) = [C(X),A(X)\, (10.5) 

( -X T {\l n -A)-^Y -X T (I n -XA)^Y\ ( I r 

[ } ~ [l r + X T (I n - \A)-iY X T (XI n - A)- X Y ) ' ^ ( ' ~ \A - Air 

A being given by J 10. J^ . 



2 For simplicity here we only give the Lagrangian description. 
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Apart from the integrals provided by the Lax matrix £(A), the equations also possess 
the matrix integral X T Y — Y T X associated to the GL(n, R)-symmetry 

(X,Y,X,Y) i — > (XR T ',YR- 1 ,XR T '.FiT 1 ), R e GL(n,R). 

Coupled Neumann system on V n , r . This systems generalizes the motion of 2 points 
x, y on the unit sphere S' 1 ^ 1 C R n that interact via the bilinear potential (x,Ay)/2. The 
latter system was introduced in [33] (see also [4"7]V 

Namely, let matrices X, Y £ M n . r (K) define two points on V n ^ r and let P,QG M njr (K) 
represent their momenta such that 

x T p + p T x = Qj Y T Q + Q T Y = 0. 
Assume that the evolution of X, Y is described by the Hamiltonian 

H = T K +tr(X T AY), (10.6) 

T K = \ tv(P T P) - Q + «) tr((jTP) 2 ) + \ tr(Q T Q) - Q + «) tr((^Q) 2 ), 

where T K is the kinetic energy of the points defined by an S'0(n)-invariant metric on V n ^ r , 
which depends on the parameter k. As above (see Section 3), for k = —1/2 we have the 
Euclidean metric and for k — the normal one. (Note that in the case r = 1 we have 
X T P = Y T Q = 0, and all the above metrics coincide.) 
The Hamilton equations with multipliers have the form 

X = P - (1 + 2k)XP t X, 
P = -AY + (1 + 2k)PX t P + ATA, 
Y = Q - (1 + 2k)YQ t Y 1 
Q = -AX + (1 + 2k)QY t Q + Qn, 

with 

A= ^(X T AY + Y T AX) -P T P, II = ^(X T AY + Y T AX) - Q T Q. (10.8) 

Borrowing the terminology of [47], we call (|10.7[) the r -coupled Neumann systems on 
V nir . They are invariant with respect to the right diagonal S l O(r)-action on the product 
T*V ntr x T*V ntr , and the corresponding matrix momentum X T P — P T X + Y T Q — Q T Y is 
preserved along their flows. 

The book [33] presented a "big" Lax pair of the r-coupled systems with the so(n 7 n)- 
matrices depending on parameter v 



(10.7) 



L{v) = [L(v),A(v)), 

XP T 

QY T - YQ T 

A\ (PX T -XP T 



Tl , fO A\ (PX T -XP T i / u .i, \ | 

m=[ A J^+( o QY T - YQ T ^ ^ x ' v 1 n )'' 

A(^) = L+(A) 






XY T 


YX T 










QY T 


- YQ T 



A OJ \ 

being zero n x n block. Curiously, this Lax pair holds only for the case of the normal 
metric (k = 0), and not for the Euclidean one, as one might expect and as happens in the 
case of the Neumann flows on T*V n:T . 

Below we also present the dual " small" Lax representation. Introduce r x r matrix 

F X (X,Y) = X T (A 2 - \ 2 )~ X Y. 
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Theorem 10.2 The coupled Neumann system \1U. 7| ) with n = 
pair with the spectral parameter X 



admits the following Lax 



±£(X) = [C(X),A(X)}, 



(10.10) 



£(A),^l(A) 



being Ar x 4r matrices 



£(A) 



/ AF A (P,X) 
-AP A (X,X) 

f a (aq,x) 



XF X {P,P) 
-XF X (X, P) 
\ r + F x {AQ,P) 
-F X (AY,P) 



F X (AP,Y) -l r + F x {AP,Q)\ 

-F X (AX,Y) -F X {AX,Q) 

XF X (Q,Y) XF X (Q,Q) 

-XF X (Y,Y) -XF X (Y,Q) J 



A(X) 



(X T P A 

I r X T P 

-AI r 

\ 



-AlA 


y t q n 

Ir Y T Q) 



and K, II are given by H0.8\) . 

Again, the proof is straightforward. For r = 1 the Lax pair (|10.10[) was given by Suris 
in @7J. 

It is still not clear whether the r-coupled Neumann system with the Euclidean metric 
admits a Lax representation. However, we can consider the following perturbation of the 
Hamiltonian (|10.6p : 



Then the corresponding flows imply the Lax representation with matrices (|10[) for any k, 
while the Lax representation (110. lip holds with X T P and Y T Q on the diagonal of -A(A) 
replaced by (1 + 2k)X t P + 2kY t Q and 2kX t P + (1 + 2k)Y t Q, respectively. In particular, 
by taking n = —1/2, we get the Lax representation of an r-coupled Neumann system with 
the Euclidean metric and with an additional interacting term tr(X T P Y T Q). 

Complex Stiefel manifolds. The Neumann systems and geodesic flows on V n , r can be 

extended to the complex Stiefel varieties W n ,r as well. Recall that W n ,r is the space of r 
ordered orthogonal vectors (z\, . . . , z n ) in C™ endowed with the standard Hermitian metric, 
or equivalently, the set of n x r matrices Z £ M„ jr (C) satisfying 



(see, e.g., [29 J. The variety W n , r can also be identified with the homogeneous space of the 
unitary group: W n , r = U(n)/U(n — r). The real Stiefel variety V n , r is thus a submanifold 
of given by the condition Z = Z . 

While the idea of integrable geodesic Sows on W n<r follows from the general construction 
given for compact homogeneous spaces [TDl [TTl [T2] , to our knowledge, potential systems on 
W n . r for r > 1 were not studied yet. 

We shall consider the Neumann system with the metric induced by the Hermitian metric 
and defined by the Lagrangian function 



H K = T K +tr(X T AY) - 2k tr(X T P Y T Q). 



■r 



(10.11) 




(10.12) 



As above, the matrix A is a real diagonal n x n matrix. 
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The Euler-Lagrange equations with multipliers read: 



Z = -AZ + ZK, Z=-AZ + Zh, (10.13) 

where 

A = Z T AZ — Z T Z = A T . (10.14) 

Then the Neumann system on V n ,r with the Euclidean metric (I5.5P can be regarded as 
a subsystem of (|10. 13|) : if Z(t) is its solution on the complex Stiefel variety W n , r with the 

initial conditions satisfying Z = Z ', Z = Z ', then (X(t),P(t)) = (Z(t),Z(t)) is a solution of 
the Neumann system (|5.5[) on V n<r , and vice versa. 
We also have 

Theorem 10.3 The Neumann system on the complex Stiefel manifold H10.13\) imply the 
2r x 2r matrix representation with the spectral parameter X 

j t £(X) = [£(X),A{X)}, (10.15) 

rn ,_(-Z T (XI n -A)- 1 Z -Z T (I n - AA)" 1 ^ An ,f I r 
L[A) \l r + zm n -XA)^Z ZT(\I n -A)-iZ A[A) -\A-XI r 

with A given by fal0.14\ l- 

The system is invariant with respect to a right U (r)-action and the symmetry with 
respect to the left action of U(l) n defined by 

Z^->djag(pi,... lPn )-Z, Pi€U(l), i = l,...,n. (10.16) 

The symmetries imply the conservation of u(n) and it(l)™ momentum maps 

r 

4>o = Z T Z - Z T Z and ipj — zfzf - zfzf, j = 1, . . . , n, 

i=i 

where Z = (zi, . . . , z r ). 

Under the condition p\ = • • • = p n , (|10.16|) dehnes the left J7(l)-action, which is free. 
Then we have a well defined reduced Neumann flow on the quotient space (complex projective 
Stiefel variety) PW n , r = W n , r /U(l). 

On the other hand, the right J7(r)-symmetry enables one to reduce the Neumann system 
to the complex Grassmann variety Gc(n,r) = U(n)/U(r) x U(n — r) of r-dimensional 
complex planes in C n , or, in general, to {/(n)-adjoint orbits O = U(n)/U(ki) x • • • x U(ki) x 
U(n — r), where ki + &2 + ■ ■ • + ki = r. 
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